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Abstract 

This work is focused on the longtime behavior of a non linear evolution 
problem describing the vibrations of an extensible elastic homogeneous beam 
resting on a viscoelastic foundation with stiffness k > and positive damping 
constant. Buckling of solutions occurs as the axial load exceeds the first 
critical value, (3 C , which turns out to increase piecewise-linearly with k. Under 
hinged boundary conditions and for a general axial load P, the existence of 
a global attractor, along with its characterization, is proved by exploiting a 
previous result on the extensible viscoelastic beam. As P < /3 C , the stability 
of the straight position is shown for all values of k. But, unlike the case 
with null stiffness, the exponential decay of the related energy is proved if 
P < j3(k), where /3(k) < f3 c (k) and the equality holds only for small values 
of k. 
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1. Introduction 



1.1. The model equation 

In this paper we investigate the longtime behavior of the following evo- 
lution problem 



d tt u + d xxxx u - f/3 + J \d(u{i, -)\ 2 d^j d xx u = -ku - 5d t u + /, 

u(0,t) = u(l,t) = d xx u(0,t) = d xx u(l,t) = 0, 
u(x, 0) = uq(x), d t u(x, 0) = ui(x), 



in the unknown variable u = u(x,t) : [0, 1] x R + — > R, R + = [0, oo), which 
represents the vertical deflection of the beam. For every x G [0, 1], uq, u\ are 
assigned data. The real function / = f(x) is the (given) lateral static load 
distribution and —ku — 5d t u represents the (uniform) lateral action effected 
by the medium surrounding the beam. Finally, the parameter e R accounts 
for the axial force acting in the reference configuration: /3 > when the beam 
is stretched, (3 < when the beam is compressed. Usually, the axial load is 
referred as P = — (5. 

The solutions to problem (II. ip describe the transversal vibrations (in di- 
mensionless variables) of an extensible elastic beam, which is assumed to have 
hinged ends and to rest on a viscoelastic foundation with stiffness k > and 
damping constant 5 > 0. The geometric nonlinearity which is involved ac- 
counts for midplane stretching due to the elongation of the bar. A simplified 
version of this beam model has been adopted to study the vibration of rail- 
way track structures resting on a viscoelastic soil (see [lj]). There, the elastic 
and damping properties of the rail bed are accounted for by continuously 
distributed or closely spaced spring-damper units. 

In recent years, an increasing attention was payed on the analysis of 
vibrations and post-buckling dynamics of nonlinear beam models, especially 
in connection with industrial applications. For a detailed overview, we refer 
the reader to [l9| and references therein. Nowadays the study of this subject 
has become of particular relevance in the analysis of micromachined beams 
|25| and microbridges [16( 



It is worth noting that the static counterpart of problem (11.11) reduces to 



u"" 



^ + J o l \d^-)\ 2 d^ u" + ku = f, 



(1.2) 

u (0) = u (l) = u"(0) = = 0. 
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Obviously, these steady-state equations does not change in connection with 
dynamical models accounting for any kind of additional damping, due to 
structural and/or external mechanical dissipation. When k = the inves- 
tigation of the solutions to (11 .2p and their stability, in dependence on (3, 
represents a classical nonlinear buckling problem in the structural mechan- 
ics literature (see, for instance, fj, 0, [io|]) which traces back to the pioneer 
paper by Woinowsky-Krieger [24|. Numerical solutions for this problem are 
available in the literature (see, for instance, (if). Recently, a careful analysis 
of the corresponding buckled stationary states and their stability properties 
was performed in 18] for all values of f3. In [3| this analysis was improved to 
include a more general nonlinear term and a source / with a general shape. 
As far as we know, no similar analysis of (II. 2p for the case k > is present 
in the literature. 

Neglecting the stiffness of the surrounding medium, k = 0, the global 
dynamics of solutions to problem (II .ip has been first tackled by Hale 14j . 
who proved the existence of a global attr actor for a general (3, relying on 
the existence of a suitable Lyapunov functional. The corresponding problem 
for an extensible viscoelastic beam has been addressed in Nl, when 5 > 0, 



and in [13J, when 5 = 0. In spite of the difficulty which is represented by 
the geometric nonlinearity, in all these papers the existence of the global 
attractor, along with its optimal regularity, is obtained for a general value 
of f3 by means of an abstract operator setting. The analysis of the bending 
motion for an extensible thermoelastic beam is even more tangled since the 
dissipation is entirely contributed by the heat equation, where the Fourier 
heat conduction law is assumed. Neverthelesss, the existence of a regular 
global ^attractor can be shown even including into (II. ip the rotatory inertia 
term 
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A common feature of previously quoted results is the following: if / = 
the exponential decay of the energy is provided when f3 > — /3 C , so that the 
unique null solution is exponentially stable if P < (3 C . On the contrary, as 
the axial load P exceeds (3 C the straight position loses stability and the beam 
buckles. So, when [3 < —f3 c a finite number of buckled solutions occurs and 
the global (exponential) attractor coincides with the unstable trajectories 
connecting them. For a general time-independent source term /, the number 
of buckled solutions may be infinite and the attractor coincides with the 
unstable set of the stationary points. The positive critical value (3 C is named 
Euler buckling load and, in the purely mechanical case, it is equal to the 
square root of the first eingenvalue of the d xxxx operator (which is referred as 
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Ai in the sequel). In the thermoelastic case, because of the thermal expansion, 
the mean axial temperature of the beam also affects the value of j3 c (see 121]). 



1.2. Outline of the paper 

At a first sight, problem (11. ip with k > looks like a slight modification 
of previously scrutinized models where k vanishes. This is partially true. In 
particular, we remark that the restoring elastic force, acting on each point of 
the beam, opposes the buckling phenomenon. So, the Euler buckling limit f3 c 
is no longer equal to a/Ai, but now turns into an increasing piecewise-linear 
function of k. When the lateral load / vanishes, the null solution is unique 
provided that (3 > —(3 c (k), and buckles when f3 exceeds this critical value. In 
general, as well as in the case k = 0, the set of buckled solutions is finite, but 
for some special positive values of k, called resonant values, infinitely many 
solutions may occur (see Theorem [1]). 

By paralleling the results for k = 0, the null solution is expected to 
be exponentially stable, when it is unique. Quite surprisingly, it is not so. 
For large values of k, the energy decays with a sub-exponential rate when 
—j3 > j3 > —j3 c (see Theorem [6]). In particular, for any fixed k > X±, the 
positive limiting value (3(k) is smaller than the critical value /3 c (h), and the 
former overlaps the latter only if < k < Ai. A picture of these functions, 
as k runs the positive axis, is given in Fig. [2j 

The plan of the paper is as follows. In Section 2 we discuss the general 
functional framework of (11. ip and exact solutions for the stationary post- 
buckling problem is presented for all k > 0, when / = 0. After formulating 
an abstract version of the dynamical problem, the existence of an absorbing 
set is addressed in Section 3. Some preliminary estimates and the exponen- 
tial stability result are established in Section 4. The main result concerning 
the existence of a regular global attractor is stated in Section 5, where the 
asymptotic smoothing property of the semigroup generated by the abstract 
problem is proved via a suitable decomposition first proposed in [13 . 



Although we assume here that both ends of the beam are hinged, different 
boundary conditions for u are also physically significant, such as when both 
ends are clamped, or one end hinged and the other one clamped. On the 
contrary, the so-called cantilever boundary condition (one end clamped and 
the other one free) is not covered because it is pointedly inconsistent with 
the extensibility assumption of the model. Nevertheless, the hinged case we 
consider here is very special. Indeed, other boundary conditions lead to a 
completely different analysis that must take into consideration very special 
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estimates for the complementary traces on the boundary, and only weaker 
forms of the regularity properties of solutions remain valid (see, for instance, 

0)- 

It is worth noting that several papers (see, for instance, 0, @, are de- 
voted to approximations as well as numerical simulations in the modelling of 
the deformations of extensible beams on elastic supports. In this connection, 
our paper which exhibits exact solutions is of interest in order to fit computer 
applications. In particular, our model can be useful for engineering applica- 
tions involving simply supported bridges subjected to moving vertical loads. 



For instance, it may be adapted to the study of lively footbridges [22|, [23 



Finally, we remark that our analysis is carried over an abstract version 
of the original problem which is independent of the space dimension, so that 
it could be extended to scrutinize shear deformations in plate models. The 
techniques of this work apply to plate models as well, without substantial 
changes. 

In addition, our strategy can be generalized to the investigation of non- 
linear dissipative models which describe the vibrations of extensible viscoelas- 
tic beams where the dissipative term derives from the internal viscoelastic 
dissipation (memory). Moreover, we are going to scrutinize the longterm 
damped dynamics of extensible elastic bridges suspended by flexible and 
elastic cables. In this model, the term —ku is replaced by —ku + and it rep- 
resents a restoring force due to the cables, which is different from zero only 
when they are being stretched. 



2. Stationary solutions 

Our aim is to analyze the multiplicity of solutions to the boundary value 
problem (11. 2p . Letting L 2 (0, 1) the Hilbert space of square summable func- 
tions on (0, 1), the domain of the differential operator d xxxx appearing in 

QE2D is 

V(d xxxx ) = {w G # 4 (0, 1) : w(0) = w(l) = w"(0) = w"(l) = 0}. 

This operator is strictly positive selfadjoint with compact inverse, and its 
discrete spectrum is given by A n = n 4 7r 4 , n E N. Thus, Ai = 7r 4 is the 
smallest eigenvalue. Besides, the following relation holds true 

(d xxxx ) 1/2 = -d xx , V{-d xx ) = # 2 (0, 1) n H%(0, 1). (2.1) 
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For every k > 0, let 

k 

H n (k) = + n 2 7r 2 , /3 c (k)=mm.fi n (k). 
Assuming that % 6 N be such that fi nk = min/i n (fc), then it satisfies 

neN 

k 

(n k -l) 2 n 2 k < <n 2 k (n k + l) 2 . 

As a consequence, (3 c (k) is a piecewise-linear function of k (see Fig. [2]). The 

set 

K = { 1 2 J V:i )J 6N,kj} 

is referred as the resonant set: when fc £ K there exists at least a value 
/Xj(/c) which is not simple. Indeed, /ij = /ij, z 7^ j, provided that k = i 2 j 2 7T A 
(resonant values). In the sequel, let fi m (k) be the smallest value of {/x n } nS N 
which is not simple. Of course, the /i n (fc) are all simple and increasingly 
ordered with n whenever k < Att 4 . Given k > 0, for later convenience let n+ 
be the integer-valued function given by 

n*O0) = Wp\, Up = {n e N : + fx n {k) < 0}, 

where \J\f \ stands for the cardinality of the set J\f. 

In the homogeneous case, we are able to establish the exact number of 
stationary solutions of ( II. 2p and their explicit form. In particular, we will 
show that there is always at least one solution, and at most a finite number 
of solutions, whenever the values of /j, n (k) not exceeding —/3 are simple. 

Theorem 1. If (3 > — /3 c (k), then for every k > system Al.fy) with f = 
has the null solution, corresponding to the straight equilibrium position. 
Otherwise: 

• if k G TZ and f3 < —fi m (k), the smallest non simple eigenvalue, there 
are infinitely many solutions; 

• whether k e TZ and —fi m (k) < < —/3 c (k), or k ^ 1Z and /3 < —(3 c (k), 
then besides the null solution there are also 2n+((3) buckled solutions, 
namely 

Un(x) = sm(nnx) , n = 1, 2, . . . , n* (2.2) 

with 

At = ±— V-2 [P + »n{k)\. (2.3) 
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Proof. Clearly, u = is a solution to (11. 2p in the homogeneous case for all 

1 

k and (3. To find a nontrivial solution u, we put h = (3 + J \u'(^)\ d£ , so 

o 

that u solves the differential equation 

u"" - hu" + ku = , heR, k>0. 
Letting A 2 = x, the characteristic equation 

A 4 - h A 2 + k = 

admits solutions in the form 

Xi,2 = ^ • ( 2 ' 4 ) 

As a consequence, taking into account the hinged boundary conditions, we 
obtain 

• if h > 2\/k, then xi,2 £ K + and all corresponding values of A are real, 
hence u = 0; 



if h < -2Vk, then xi,2 G K~ and Ai, 2 = ±y/\xT\i, A 3i4 = iylxd 
hence 

u = 2 i a sin Wi x + 2 2 6 sin 0J2 x, 



where a and b are suitable constants, while 



y^^nvr, W2 = v ^ = £7r, (2.5) 



Now, letting 2 a = — %A and 2 6 = — ?5 , we can write the solution into the 
form 

u = A sin nir x + B sin £ir x. (2-6) 



Moreover, from ( 12. 4ft . ( 12. 5 j) we obtain 



1 

^(-h-Vh 2 -4k) = nV /i = -Hn(k) = (3 + J |w'(0| 2 ^ (2.7) 
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1 

I ( - h + y/h 2 -Ak) = fir 2 =>h = -y. t {k) = (3 + f | 2 d£ (2.8) 



o 



from which it follows fie(k) = fi n (k). In order to represent J °^ m 

o 

explicit form, we are lead to consider two occurrences. 

• Let k 1Z. In this case, from (12 .7p and (I2.8P it follows £ = n and, by 
( 12. 6p . the following equality holds 

i 

/ 1^(0 1 2 d£ = \ [A + B) 2 nV = -f3 - //„(&), 
o 

so that, letting A = A + B, (12. 3p can be easily obtained. Of course, 
such nontrivial solutions exist if and only if —(3 > f3Jk) = min //„(&). 

neN 

• Let k G 1Z. Then, fie(k) = fi n (k) for some £ < n. In this case, from 
(12.61) we obtain the equality 

i 

/ |u'(£)| 2 d£=iiW + ~£W, £^n, 
o 

and (I2.7l) - (l2.8p cannot uniquely determine the values of A and B. 
Accordingly, (12. 6p represent infinitely many solutions provided that 
-p > n m {k). 

Remark 1. Assuming k = we recover the results of 18] and Q. 

When k (jL 1Z the set of all stationary states is finite and will be denoted 
by Sq. Depending on the values of k and (3, the pairs of solutions branch 
from the unbuckled state A^ = at the critical value (3 = —f3 c (k), i.e., the 
beam can buckle in either the positive or negative directions of the transverse 
displacement. These branches exist for all < —/3 c (k) and are monotone 
increasing functions of \(3\. For each n, (I2.3P admits real (buckled) solutions 
A^ if and only if (3 < —fi n . When k < 47r 4 , for any (3 in the interval 

h h 

™ ' - \1 2 ^ a , 22 



(n + lY Tl 2 < /3 < =—r - ?T 7T 



(to + l) 2 7r 2 to 2 7r 
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.4 




the set Sq contains exactly 2n* + 1 stationary points: the null solution and 
the pairs of solutions represented by (12.21) . These properties are sketched in 



FigH](see also |18|). 



3. The Absorbing Set 

We will consider an abstract version of problem (II. ip . To this aim, let H be 
a real Hilbert space, whose inner product and norm are denoted by (-, ■) and 
|| • ||, respectively. Let A : T>(A) d H — > H be a strictly positive selfadjoint 
operator. We denote by Ai > the first eigenvalue of A. For £ G K, we 
introduce the scale of Hilbert spaces 

H t = V(A e / A ), (u,v) t = {A"l\A^v) } \\u\\, = \\A e /\\\. 

In particular, Hi + \ d Hi and the following scale of Poincare inequalities 
holds 

>Ai INI/ < ll^lll+i- (3-1) 
Finally, we define the product Hilbert spaces 

He = Hg +2 x Hi. 
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The abstract problem 

For G K. and / G Hq, we investigate the following evolution equation 
on T-Lq 

d tt u + Au+(f3+ \\u\\ 2 l )A ll2 u + Sd t u + ku = f, (3.2) 

in the unknowns u(t) : [0, oo) — > H2 and dtu(t) : [0, 00) — > Hq, with initial 
conditions 

(u(0),d t u(0)) = (u ,ux) = z G "H - 

Problem (II. ip is just a particular case of the abstract system (I3.2p . obtained 
by setting H Q = L 2 (0, 1) and A the realization of d xxxx in H . 

Eq. (13. 2 p generates a strongly continuous semigroup (or dynamical sys- 
tem) S(t) on 7-L Q : for any initial data z G Ho, S'^Zo is the unique weak 
solution to (13. 2p . with related (twice the) energy given by 

£(t) = \\S(t)z \\ 2 Ho ^ \\u(t)\\l+\\d t u(t)\\ 2 - 

Besides, S(t) continuously depends on the initial data. We omit the proof 
of these facts, which can be demonstrated by means of a Galerkin procedure 
(e.g., following the lines of Q). The crucial point in applying this technique 
is to have uniform energy estimates on any finite time-interval. As it will be 
apparent, these estimates are easily implied by the uniform inequalities on 
(0, +00) established in the subsequent sessions. 

Now, we prove the existence of the so called absorbing set for the flow 
generated by problem (13. 2p . that is, a bounded set into which every orbits 
eventually enters. Such a set is defined as follows: 

Definition 1. Let B(0, R) be the open ball with center and radius R > 
in "H - A bounded set B^ C "Ho is called an absorbing set for the semigroup 
S(t) if, for any R > and any initial value Zq G -8(0, R), there exists t (R) > 
such that 

z(t) G B no V t > t , 
where z{t) = S(t)z is the solution starting from zq. 

The main result of this section follows from two lemmas, involving the 
functionals £(t) and $(i) defined as 

C(t) = £{t) + ^+\\u{t)\\\) 2 + k\\u{t)\\ 2 >£{t)>^ (3.3) 

$(t) = C(t) +e(d t u, u). (3.4) 
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Lemma 2. For all t > and z G "Ho with ||zo||"H — ^> i/iere exzsfo a 
positive constant C (depending on \\f\\ and R) such that 

€(t) < C. (3.5) 

Proof. If we consider the functional 

T(t) = £(*)- 2 </,«(*)), 

from the energy identity 

_ = _2^(/3 + ||u||J)^5u + jbw + (jg t u-/, d t u 
it easily follows the decreasing monotonicity of T 

7 -7— 7 

— = — + 2 (0 + \\u\^) (a\u , A^d t u) + 2k (u, d t u) -2 {f, d t u) 
= -28\\d t u\\ 2 <U 

and then 

m<m<c (R, n/ii). 

Taking into account that 

IKOII 2 <^IKOIl5<^(*) = Cif(0, 



we obtain the estimate 



.F(f) > £{t)-2 (f , «(*)> > Sit)- 1 - \\f\\ 2 -e \\u(t)f > (l-eCO^f)-- 
Finally, fixing e < ^, we have 

< tztct^ ) + \ ii/ii 2 ) ^ r^W (c, ° (i? ' + \ ll/l|2) = a 

Lemma 3. For any given z G "Ho orad /o r i > and (3 G R, w/ien 5 is 
sraali enough there exist three positive constants, mo, m± andm 2 , independent 
oft such that 

m £(t) < $(t) < m 1 £(t) + m 2 . (3.6) 
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Proof. In order to prove the lower inequality we must observe that, by 
Young inequality 



I {dtu ,u)\>- 



1 



2 1 II 1 1 2 

"2 H-ll 



hence, we obtain 

m > \\u(t)\\i+(i-~) \\d t u(t)\\ 2 +~(p+\\u(t)\Q 2 +(k-~) Mm 2 . 

If we choose e small enough to satisfy e < 2 and e < 2k, then we have 

$(t) >m £(t) >m S(t), (3.7) 

where m = min{l — |, 1 — ^}. 

The upper inequality can be obtained using the definition of $ and ap- 
plying the estimate 



\(d t u, u)\ < 



1 



2 , 1 ll n 2 

+ 2A-NI 2 



First, we can write 
<&(*)< 



(3.8) 



\Ht)\\l+(^ + i)\\dMt)\\ 2 + ^ + \Ht)\\l) 2 ■ 



Then, by (13. 1 p and Lemma [2] we infer 



{P + \Wl)<\p\ + 



c = c, 



(3.9) 



so that we finally obtain 
*(*)< 



+ 2^ = m i^W + m 2 • 



Theorem 4. For any /5 6 K, there exists an absorbing set in Ho for the 
dynamical system (S(t), / Ho). 

Proof. By virtue of (I3T3"|) we have 

^ = ^ + 2 (/3 + \\u\Q (j&u , + 2k (u, d t u) = 



-25 



2{f,d t u). 
12 



Moreover, by (13. 4p 

^ = ^ + e (u , d tt u) + e ||^f = - (25 - e) \\d t u\\ 2 + 2(f, d t u) + 
at at 

+ e [— \\u\\1 — 5 (d t u, u) — k \\u\\ 2 — (3 \\u\\1 — \\u\\* + (f , u 
A straightforward computation leads to the identity 



^+6$ + 2(5-e) \\d t uf + ~\\u\\\ = 

= 2 (/ , d t u) +e[(f,u)-(5- e) (d t u, it)] + |/3 2 . 



(3.10) 



In the following we estimate the terms in the rhs of the previous equality. 
By means of the Holder and Young's inequalities and (13. 8p . we have 

d ^ + e$ + ±(35-5e)\\d t u\\ 2 < 



Now, choosing e < I 5 and e < 1 + 5 we find 



( ^--s<t><ls3 2 -(--l) f > 



dt - 2 ^ V £ 2 

where 

m = — (l + £-e) >0. 

Finally, from Lemma [3] we obtain 

1 2 /l 1 
dt ~ 2 ' \e 2 

where w — e ( 1 — — £ ) is positive provided that e < — . Since $ is positive 

\ mo / ?n r 

also when (3 < provided that £ is chosen small enough, using Gronwall 
lemma it follows that 

$(*)<$ (0) e — ' + i£/3 2 + Q + HZ || 2 , (3.11) 
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and accordingly 

E {t) < J-$(f) < r (i?)e- rot + rxCfl, 11/11), (3.12) 

where 

^ + (i + i)n/ir 

As a consequence, every ball £> (0, R) in "H with radius R > 1 + Ti(f3, \\f\\) 
can be chosen as an absorbing set in that it verifies the following statement: 
for all Zq = (u , Mi) G B(0,R), there exists t (R) = ^logr (i?) such that 
for any t > t , z(t) G B(0, R). 



T (R) 



m 



and ^08,11/11) = — 



1 

m 



4. Exponential Stability 

A direct proof of the exponential decay of the energy seems out of reach, 
so we exploit the equivalence between the energy 8 and the functional 

which can be proved to be exponentially stable. The positivity of such a 
functional will be obtained as a direct corollary of the next Lemma EJ 

Recalling Theorem [TJ the set S of stationary solutions reduces to a sin- 
gleton when 



(3 > -f3 c (k) = - min jj n (k), fi n (k) = V A n 



ngN 



1 + f 



A n = n 4 vr 4 . (4.1) 



It is worth noting that (3 c (k) is a piecewise-linear function of k, in that 
(3 c {k) = [ii{k) when < k < V%\/A2 = 47r 4 , and in general 



c (k) = n n (k) when a/A„_ia/A^ < k < V\A+i' 

Unlike the case k = 0, the energy £{t) does not decay exponentially in the 
whole domain of the (/3, k) plain where (14. ip is satisfied, but in a region which 
is strictly included in it. 
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For futher purposes, let 

f3 c (k) < k < Ai, 

2Vk k > Ai. 

A picture of this function is given in Fig. |2j 
Lemma 5. Let (3 G K., > and 

Z/u = Am + + £;« . 

There exists a real function v = z/(/3, k) such that 

(Lu , u) > v 1 1 1 1 g , 
where v((3, k) > if and only if (3 > —/3(k). 
Proof. Taking into account the inner product 



we put 
so that 



(Lu , u) = IHI2 + /^IMIi + ^IMI 2 



X = \\u\\ 2 ; Y = \\u\li ; Z = \\u\ 



(Lu , u) = I(X, Y,Z)= X 2 + /3 Y 2 +k Z 2 . 
If (3 > 0, the desired inequality trivially holds true by choosing v 

1 1 2 

Letting p < 0, by means of the interpolation inequality Wu^ < 
we obtain 

I(X, Y, Z) > X 2 + (3XZ + kZ 2 = J(X, Z). 
Hence the thesis can be rewritten as follows: find v > such that 

J(X,Z) \ Do >vX 2 , 

where D = { (X, Z) : X > 0, Z > 0, < Z < X/^/\[} . 

In order to prove this statement, we introduce the set 

M = <!m Gl: 0<m< 1 
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so that D = {(X, Z) : Z — ml, X > and m G M} and the original prob- 
lem reads: find v > such that 



J(X, Z) | z=rnX > v m X 2 , Vm G M 



and 



v = inf v m > 0. 

ra€M 



We first observe that 



J(X, Z) \z=mX 



(1 + l3m + km 2 )X 2 = rj{m)X 2 



where 77(771) is a concave parabola. Hence, we have to find the region in the 
(/3, fc)-plane where 77 admits a strictly positive minimum on M, that is to say 
Vmin = v > 0. We shall prove that v > if and only if /3 > —(3(k). To this 
end, we split the proof into four steps. 



• Step 1. We consider the region 

#1 = |(/3, k) : k > 0, -2Vk < (5 < j. 

Since the discriminant of the parabola is negative, = 1 — 4k < 0, 
the value 77(771) is strictly positive for all m in the closed interval M, 
and then u > 0. 



• Step 2. Let i? 2 = <^ (/3, fc) : < k < Ai, -v 7 ^ - k/y/\[ < (5 < -2^k 



Observing that Ai < ^/Xi^/Xi, we infer that f3 c (k) = fii = \f\[+k/ \f\[. 
Now, since A v = j3 2 — Ak > 0, there exist two solutions, mi,m 2 G R, 
of 77(777) = so that i] changes sign on R. In particular, it must be 
negative inside the open interval (mi , 7772) because of k > 0. On the 
other hand, 77 is positive at the ends of M. Indeed, 77(0) = 1. In order 



to evaluate the sign of 7](l/v%") we let k — p\i, with < p < 1, and 

we obtain 





According to the definition of R2 




(4.2) 



therefore 




+ (l+p) >0, 
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which implies 77 fl/\/Ai) > 0. Thus we infer that either (mi , m 2 ) C 
M, or (mi , 772,2) is external to M. In the sequel we prove the latter 
occurrence by showing that the vertex of the parabola lies outside of M. 
Letting m* the abscissa of the vertex, it satisfies 7/ \ m=rrlf = f3 + 2km* = 
so that 





772* 



2k 2pAi 

So, in order to have 77(777) \ M > 0, it is enough to prove that 

m* > 



2pAi v 7 ^ 
which is trivially true because from (14. 2ft we have 

Step 3. Let # 3 = |(/3 , ib) : A: > Ai, -\Ai~ < P < -2^}. 

When (3 < —2\/k, as before we have A v = (3 2 — Ak > 0, so that 77 is 
negative valued in the open interval (mi , 7772) delimited by solutions 
mi , r?72 to the equation 77(777) = 0. Nevertheless, in this case the vertex 
of the parabola lies inside M. Indeed, if we take k = a Ai, with a > 1, 
the abscissa of the vertex satisfies 

1 
< m* = — - — — < 



2aAi ' 

which holds true by virtue of the definition of -R3, in that 
-/?<(! + a) < 2o\f\ x . 



As a consequence, the minimum of 77 on M is rj (m*) < 0. When 
= —2\fk, we have mi = 7772 = 777* and 7/(777*) = 0, so that the mini- 
mum of rj on M vanishes. In both cases, ?7 m i n = v is not positive. 



• Step 4. We consider the set R 4 = { (0, k) : k > 0, < -\Ai - /c/\Ai}- 
In this case 7/(0) = 1, whilst 

"(^t)^(^ + 5 + ^7) <0 

and the minimum of 77 on M cannot be positive. 



17 



We are now in a position to prove the following 

Theorem 6. When f = 0, the solutions to U.l\) decay exponentially, i.e. 

S{t) < c S{0)e- ct 

with Co and c suitable positive constants, provided that (3 > —fl(k). 

Proof. Let $ the functional obtained from $ by letting / = and neglect- 
ing the term |/3 2 , i.e. 

<5(t) = ||«(t)||^+||^(t)|| 2 +/3 ||«(t)||;+i \\u(t)\\i+k \\u(t)\\ 2 +e (d t u(t) , u(t)> . 
In view of applying Lemma O we remark that 

$ = (Lu , u) + \\d t u\\ 2 + - \\u\\\ + e (u, d t u) . 

The first step is to prove the equivalence between £ and $, that is 

c\ £ < # < c 2 £ . 
We split the proof into two parts. 
• First step: c\ £ < $ 

By virtue of Lemma [5] and (13.11) . the following chain of inequalities holds 
provided that (3 > —j3(k), which ensures the positivity of v. 

£ (» - 4) + " I ) 119,1,112 * min {" " n? 1 " 1} £ ' 

If we choose e < min {2v\i , 2} and we put C\ = min jz/ — 1 — f |, it 
follows 

Ci £ < $. 
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• Second step: $ < c 2 £ 



Using the expression of $, Young inequality and ( 13. ip . we can write the 
following chain of inequalities 



$ < 1 



k 1 

+ 



Ai 2Ai 

< 



\u\ 



1 + 



k 1 
2 + — + 



Ai 2Ai 2 
In particular, from (I3.9P we obtain 



2 , a II II 2 , ^ || ||4 ^ 

+ P \\ u \\i + 2 W u Wi - 



1 



f- K + /?!!< +2 Nit- 



$ < 2 + 



k 



1 



Ax 2Ai 



k 1 e 2 



£ + C \\u\\{ < 
C 



£ = c 2 £. 



The last step is to prove the exponential decay of $. To this aim, remember- 
ing the hypothesis / = 0, we can write the identity (I3.10p in the following 

way: 

Exploiting the Young inequality we obtain 



+ 2 (5 - e) 



2 + | IMIi +s{S — e) (d t u , u) = 0. 



and choosing £ < 5 it follows 



< Hl 2 



2Ai 



d$ - s 2 (S-e) „ e 2 (5 - e - 

— + e$< — ^- - £ < K - ' $ 

dt ~ 2Ai _ 2ciAi 



and finally, if e is small enough, we have 

— + c$ < 
at 



(4.3) 



where c = e 



1 - 



e(g-e) 
2ciAi 



> 0. Equation ( 14. 3ft implies 



ci£(t) < * (t) < $ (0) e- ct < c 2 £ (0) e 



-ci 



The thesis follows by letting cq 



£1 

ci 
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Remark 2. We stress that Theorem holds even if k = 0. In this case how- 
ever we have (3(0) = — /3 C (0) = — y/X%. Then the null solution is exponentially 
stable, if unique. 

In Figj2j the piecewise-straight line is the bifurcation line (3 = —j3 c (k): 
below it, the system has multiple stationary solutions, above it, there exists 
only the null solution. The curve (3 = —(3(k) is composed of the straight 
segment (3 = /ii(fc), when < k < A 1; and the parabola (3 = —2\/k, if 
k > Ai- It is worth noting that each segment composing the graphic of — (3 C 
is tangent to the parabola at k = \ n = n 4 7r 4 , n e N. In the region of the 
plain (f3, k) which is bounded from below by (3 > —(3(k) the exponential 
stability holds true. Whilst, in the area between the red and the blue line, 
when k > Ai, we have asymptotic, but not exponential stability. 

e\ 



•^1 4Ai A2 




Figure 2: A picture of the functions —f3(k) and —(3 c (k). 



5. The Global Attractor 

We now state the existence of a global attractor for S(t), for any (3 EM. and 
k > 0. Recall that the global attractor A is the unique compact subset of 
H which is at the same time 
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(i) attracting: 

Hm6{S(t)B,A)->0, 

t— >OD 

for every bounded set B C %o, 

(ii) fully invariant: 

S(t)A = A, 

for every t > 0. 

Here, <5 stands for the Hausdorff semidistance in H , defined as (for Bi, B 2 C 
•Ho) 

5(i3i,i3 2 )= sup inf \\z 1 -z 2 \\ Ho . 



We address the reader to the books [14|, |2l| for a detailed presentation of the 
theory of attractors. We shall prove the following 

Theorem 7. The semigroup S(t) on Ho possesses a connected global attrac- 
tor A bounded in H2 = D(A) x D(A^). Moreover, A coincides with the 
unstable manifold of the set S of the stationary states of S(t), namely, 

z is a complete (bounded) trajectory of S(t) : 
Z{0): £m || z(-t) - S\\ Ho = 

The set S consists of all the vectors of the form (u, 0), where u is a (weak) 
solution to 

Au+ (P+\\u\\\)A l l 2 u +ku = f. 
It is then apparent that S is bounded in If S is finite, then 

z is a complete (bounded) trajectory of S{t) 
A = I z (0) : sucri that 3 zi, Zi G 5" : ^ . (5.1) 

hm - ^|| Wo = hm \\z(t) - z 2 \\ Hq = 



Remark 3. When / = 0, k > and > -p c (k), then A = S = {(0,0)}. If 
< —/3 c (k), then S = Sq may be finite or infinite, according to Theorem [TJ 
In the former case, (15.11) applies. 
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The existence of a Lyapunov functional, along with the fact that S is 
a bounded set, allows us prove Theorem [7J by paralleling some arguments 



devised in 13 



By the interpolation inequality \\u\\1 < ||u||||u||2 and (13. ip . it is clear that 

2IMI2 < Nll + /3|Hli + 7lHI 2 <"i|Hl2, (5-2) 

for some m = m(/3,7) > 1, provided that 7 > is large enough. Now, 
choosing 7 = ji + k, where k > is a fixed value, we assume /1 large enough 
so that (15.21) holds true. Then, we decompose the solution S(t)z into the 
sum (see |l3| ) 

S(t)z = L{t)z + K{t)z, 

where 

L{t)z = (v(t), d t v{t)) and K{t)z = (w(t), d t w(t)) 

solve the systems 

'd tt v + Av + (ft + WuWDA^-^v + fiv + 5d t v + kv = 0, 
(v(0),d t v(0)) = z, 



(5.3) 



(5.4) 



and 

d tt w + Aw + (j3 + WuWDA 1 ^ — fiv + Sd t w + kw = f, 
(w(0),d t w(0)) = 0. 

Having fixed the boundary set of initial data B(0,R), Lemma |2] entails 

sup {\\u(t)\\ 2 2 + \\d t u(t)\\ 2 } <C = C(R), Wz e B(0,R), 

t>o 

where (u(t) , d t u{t)) = S(t)z ; this formula will be used many times in the 
subsequent proofs. 

In order to parallel the procedure given in 13], we shall prove Theorem 
[7] as a consequence of the following five lemmas. We start showing the ex- 
ponential decay of L(t)z (see Lemmas M and [9]) by means of a dissipation 
integral (see Lemma fTUl) . Then, we prove the asymptotic smoothing prop- 
erty of K(t), for initial data bounded by R (see Lemma [TTT) . Finally, by 
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collecting all these results, we provide the desired consequence (see Lemma 

□2D- 

Henceforth, let R > be fixed and ||^||-h < R. In addition, C will denote 
a generic positive constant which depends (increasingly) only on R, unless 
otherwise specified, besides on the structural quantities of the system. 

Lemma 8. There is u = u(R) > such that 

\\L(t)z\\ Ho < Ce- Ut . 

Proof. Denoting 

S (t) = \\L(t)z\\i o = \\v(t)\\l + \\d t v(t)\\ 2 , 
for e > to be determined, we set 

%(t) = S (t) + (5\\v{t)\\l + (» + k)Mt)\\ 2 + NOIIilKOII? + e(d t v(t)At))- 



In light of Lemma [2] and inequalities (I5.2p . assuming that e is small enough, 
we have the bounds 

~£o < $o < CS . (5.5) 

Now, we compute the time-derivative of $o along the solutions to system 
( 15. 3 p and we obtain 

^ + e% + 2(5 - e)\\d t v\\ 2 = 2(d t u, A^u) \\v\\\ - e{8 - s)(d t v, v). 
at 

Using (I3.5p and assuming e small enough (in particular, e < 5), we control 
the rhs by 

C\\d t u\\£ + ~£o + (5-e)\\d t v(t)\\ 2 . 



$o + (5-^)||^|| 2 <C||^i||$ 



o- 



So, from (15. 5p . we obtain 

d®n e 
— - + - 
dt 2 

Finally, the functional $o fulfills the differential inequality 

^ + £ -^o<C\\d t u\\%- (5.6) 
The desired conclusion is entailed by applying the following two lemmas. 
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Lemma 9. (see Lemma 6.2, [131 ]) Let <p : [0, oo) — > [0, oo) satisfy 

(p' + 2ecp < gcp, 
for some e > and some positive function g such that 



J g(y)dy<c + e(t-r), Vre[0,t], 

with Co > 0. Then, there exists C\ > such that 

V?(t) < C!^(0)e" Et . 

Lemma 10. For any o > small 

t 

\\d t u(y)\\dy<a(t-r) + ^ , 

cr 

T 

/or ei>en/ t > r > 0. 
Proof. For e e (0, 1], we set 

v& = $-2(/, M )-i/3 2 , 

where $ is defined by (13 .4p . Taking the time derivative of \l/ and using (I3.10p . 
we find 

^ + e^ + £ -\\u\\t + 2(6-e) \\d t u\\ 2 = s [(/ , u) + (5 - e) (d t u , u)) . (5.7) 
By virtue of Lemma [2j ^ is bounded and hence 

^ + eV + (5-e) \\d t u\\ 2 <eC . 
at 

Since \1/ is uniformely bounded by Lemma [3j we end up with 

^ + (5-e) \\d t u\\ 2 <e[C-^]<eC. 

Integrating this inequality on [r, t] , we get 

t 

(5-e) [ \\d t u(y)\\ 2 dy<eC(t-T) + ^(r)-^(t). 
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Assuming e < |, a further application of Lemma |3] entails 

t 

j \\d t u(y)\\ 2 dy< e d(t-r) + C 2 . 

T 

Finally, by the Holder and the Young's inequalities 



t 



\\d t u(y)\\dy<Vt^(eC 1 (t-r) + C 2 y < 

< 2 ^7C[ (t - r) + = a(t - r) + ^ 

where cr = 2y/eC\. 

The next result provides the boundedness of K{t)z in a more regular 
space. 

Lemma 11. The estimate 

\m)4u> < c 

holds for every t > 0. 
Proof. We denote 

m = \\m4k = \Ht)\\i+\\d t w(t)\\i 

For £ > small to be fixed later, we set 

$i = £i + (p + \\u\\l) \\w\\l + e(d t w, Aw) - 2(f, Aw) + k\\w\\ 2 2 . 
The interpolation inequality 

\\w\\l < \\wW2WwW4 

and the fact that || if || 2 < C (which follows by comparison from (13.51) and 
Lemma [8]) entail 

/3|kll^>-^i-c. 

25 



Therefore, provided that e is small enough, we have the bounds 

~£ x - C < $! < CE X + C. 
Taking the time-derivative of $1, we find 

= 2<d t u, A 1,2 u) \\w\\l + 2fi(A 1/2 v, A 1/2 d t w) + 
+ e \fi(A 1/2 v, A 1/2 w) -(S-e) (d t w, Aw) - (/, Aw) 

Using ( I3.5p and the above interpolation inequality, if e is small enough, we 
control the rhs by 

§£ a + C\fE x + C < -E x + - < + -. 

8 4 £ 2 e 

Hence, if e is fixed small enough (in particular e < 5), we obtain 

Since $i(0) = 0, from the Gronwall lemma and the controls satisfied by $i, 
we obtain the desired estimate for E\. 

By collecting previous results, the following Lemma can be applied to 
obtain the existence of the attractor A. 



Lemma 12. (see [13], Lemma 4.3) Assume that, for every R > 0, there exist 
a positive function ipR vanishing at infinity and a compact set KLr C Ho such 
that the semigroup S(t) can be split into the sum L(t) + Kit), where the 
one-parameter operators L(t) and K(t) fulfill 

||^(*)^|ko < ^-r(*) and G JCr, 

whenever \\z\\u < -R and t > 0. Then, S(t) possesses a connected global 
attractor ^4, which consists of the unstable manifold of the set S. 
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